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ABSTRACT 

We used the Extended Submillimeter Array (eSMA) in its most extended 
configuration to investigate the innermost (within a radius of ~ 290 i?* from the 
star) circumstellar envelope (CSE) of IRC+10216. We imaged the CSE using 
HON and other molecular lines with a beam size of 0."22 x 0."46, deeply into 
the very inner edge (~ 15 i?*) of the envelope where the expansion velocity is 
only ~ 3 km s~^. The excitation mechanism of hot HCN and KCl maser lines 
is discussed. HCN maser components are spatially resolved for the first time on 
an astronomical object. We identified two discrete regions in the envelope: a 
region with a radius of < 15 i?*, where molecular species have just formed and 
the gas has begun to be accelerated (region I) and a shell region (region II) with a 
radius of 23 i?* and a thickness of 15 i?*, whose expansion velocity has reached 
up to 13 km s^^, nearly the terminal velocity of 15 km s~^. The Si^^S line 
detected in region I shows a large expansion velocity of 16 km s~^ due to strong 
wing components, indicating that the emission may arise from a shock region 
in the innermost envelope. In region II, the P. A. of the most copious mass loss 
direction was found to be ~ 120 ± 10°, which may correspond to the equatorial 
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direction of the star. Region II contains a torus-like feature. These two regions 
may have emerged due to significant differences in the size distributions of the 
dust particles in the two regions. 

Subject headings: circumstellar matter - masers - stars: AGB and post-AGB - 
stars: mass loss - stars: winds, outflows - individual(IRC+10216) 



1. INTRODUCTION 

IRC+10216 (Neugebauer &; Leighton 1969) is the most studied carbon star. Because of 
proximity (distance of 135 pc (Le Bertre 1997)), it shows a very rich spectrum of molecular 
and atomic lines particularly at millimeter and submillimeter wavelengths (e.g., Groesbeck 
et al. 1994). As such, IRC+10216 is an important target for circumstellar and interstellar 
chemistry. The carbon star is a red giant that is evolving toward the post asymptotic giant 
branch (AGB) phase. It has a high mass- loss rate of up to 10^^ M0yr~^, and is therefore 
enshrouded by optically thick circumstellar shells. 

Observing maser transitions of various molecules towards evolved stars is a powerful tool 
to study their innermost circumstellar envelopes, in particular the recent mass-loss activities 
and history of such stars. Silicon monoxide (SiO), water (II2O) and hydroxyl (OH) maser 



of oxygen rich evolved stars ( 


Morris & KnaPD ] 


976; 


Reid & Muhleman 


1978; 


Moran et al. 


1979; 


Diamond & Kemball 


2003; 


Shinnaea et al 


20041). 



Hydrogen cyanide (HCN) is one of the most abun dant molecular species in the circum- 
stellar envelopes of carbon stars (IRidgway et al.lll978l ) and is known to show maser action 
in various vibrational states (e.g., Ziurys & Turner 1986 for f > 1; Izumiura et al. 1995 for 
v=0). Since HCN is a linear triatomic molecule, it has three fundamental vibrational modes, 
(vi, V2, V3); Vi is the C-H stretching mode, V2 is the doubly degenerate bending mode, and 
V3 is the C-N stretching mode. Many observational studies have been done using HCN 
transitions in the circumstellar envelope of IRC+10216 and other carbon stars. Amazingly, 
some of the HCN maser lines are in very high energy levels (4000 K or above), far above 
than the temperature of the photosphere or the dust temperature of the region. The HCN 
component in the ^=(0, 1^^, 1) state, which is one of the high energy lev el transitions, is 
known to perturb w ith the ^=(0, 4°, 0) state through Coriolis interaction (IHocker fc Javan 
19671 ; IChantrylll97ll ). first empirically identified by Gebbie et al. (1964). The (0, 1^^, 1) and 
(0, 4°, 0) system have a high amplification in the 336.5 and 310.9 /im lines due to the fac t 
that the vibrational transition moments are sharply peaked at J = 10 (ILide fc Makil 119671 ). 
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These rovibrational transitions were detected towards IRC+10216 (jSchilke fc Menterul2003l ). 



This paper reports detections of two other hnes at such high energy levels. 



Lucas fc Guillotead (Il992l ) used the Plateau de Bure interferometer to study the details 
of the inner circumstellar envelope of IRC+10216 with the v={0, 2°, 0) J = 2 — 1 transition 
at 3 mm and they confirmed that the emission may be weak maser emission occurring in 
an acceleration region of the circumstellar envelope. However, the components of the maser 
were not significantly resolved with their beam size of 2.4"x 1.1". 

The Extended Submillimeter Array (eSMA) is a heterogeneous array that combines 
the signals detected with the eight 6-meter Submillimeter Array (SMA) antennas, the Cal- 
tech Submillimeter Observatory (CSO) 10.4 meter telescope, and the James Clerk Maxwell 
Telescope (JCMT) 15 meter telescope, operating at short millimeter and submillimeter wave- 
lengths. Including the CSO in the interferometer, one can extend the longest baseline up 
to ~782 meters, which improves the spatial resolution by ~40 %. The eSMA is one of the 
instruments that allows us to study the gas distribution and the kinematics of the inner- 
most circumstellar envelope of this carbon star in detail where the acceleration of the gas 
is taking place. Here we report an observational study on the properties of the innermost 
circumstellar envelope of IRC+10216 using HCN and other molecular lines with the eSMA. 
We detected a few hot HCN lines including the one in the ^=(0, 1^^, 1) state mentioned in 
previous paragraph. 



2. OBSERVATION AND DATA REDUCTION 

We observed IRC+10216 with HCN J = 3 — 2 transitions in various vibrational states 



along with other molecular transitions with the eSMA (iBottinelli et al.ll2008l ) on 2008 April 
14 (UT). The total bandwidth of ~ 4 GHz (2 GHz from each sideband, i.e., 264.357 - 
266.337 GHz and 274.357 — 276.338 GHz) was covered with a frequency resolution of 812.5 
kHz0. The data were taken under moderate weather conditions with zenith opacities of r ~ 
0.15 - 0.18 measured at 225 GHz. Atmospheric opacities at the observing frequency were not 
measured during the observation, however, they are close to what we measured at 225 GHz 
because these frequencies are close. Integration time for one data point was 30 seconds. The 
track lasted 5.5 hours, and about 4 hours of on-source data were acquired. The SMA was 
configured in its very extended array (maximum baselines to 508 m); the eSMA then covered 
baselines ranging from 25 to 782 m (~ 20 to 692 kA or ~ 0."2 to 10"). The phase tracking 



^Note that JCMT's receiver had a more hmited frequency coverage from 264.770 to 266.338 GHz and 
from 274.357 to 275.927 GHz. 
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center was set at RA = 9'^ 47°^ 57.^390, Decl. = 13°16'43."898 (J2000). The synthesized 
beam size was 0."46 x 0."22 with a position angle of 39.9°. 

Data reduction and cahbration were done using the MIR packag^ including passband 
and gain calibrations as well as continuum subtraction. IRC+10216 was observed together 
with quasars 3C273 and 1058+015. Passband calibration was done using the afore mentioned 
quasars. 

To make the eSMA observations, we used half wave plates (HWPs) at the JCMT and 
at the CSO to try to match the polarization angles of the receiver mixers with those at the 
SMA antennas. Based on an aperture efficiency of 75% of the SMA antennas, the estimated 
aperture efficiencies of the JCMT and the CSO telescopes were found to be 40% and 30%, 
respectively. The CSO's aperture efficiency is typically > 60% when the telescope is used as 
a single-dish telescope. The degraded efficiency of the CSO telescope was due to a failure 
of the computer that controlled the HWP. The polarization angle of the CSO's mixer on 
the sky was not matched to the SMA's. The receiver temperatures of the SMA antennas 
and the CSO antenna ranged from ~ 85 to 140 K, while the JCMT antenna's receiver 
temperature was ~ 290 K. The double sideband system temperature of the SMA antennas 
and the CSO telescope were ranging between ~ 200 - 500 K, while the double sideband 
system temperature of the JCMT telescope was ranging between ~ 300 - 600 K. 

While the continuum emission of IRC+10216 is marginally resolved on the longest base- 
lines (see Section IXT]) . we have performed self-calibration of the visibility phase gains using 
the continuum, approximated by a point source, as our primary phase calibration step. Flux 
calibration was achieved by setting the continuum component of IRC+10216 to 1.18 Jy, con- 
sistent with well calibrated observations obtained in the SMA compact array in March 2008 
whose beam size was about 3". Based on the estimated flux of 1.18 Jy for IRC+10216, the 
calculated flux of the quasar 1058+015 was 2.53 Jy, which was 9 % higher than the flux in the 
calibrator catalog measured on the same date (2008 April 17). Since the quasar is linearly 
polarized by ~ 25 - 30 % at this frequency, the 9% difference may be due to polarization. 
Note that, by applying self-calibration using the continuum component of IRC+10216, the 
data taken on the quasar 1058+015 showed that atmospheric decorrelation was eliminated 
signiflcantly even on long baselines. 

Image processing including clean was done using Astronomical Image Processing System 



^http://cfa- www. harvard. edu/^cqi/mircook. html 

■^http:/ /smal.snia.hawaii.edu/caUist/caUist.html?data=1058%2B015 
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(AIPS0) with natural weighting of the visibihties. Natural weighting was chosen to achieve 
a high signal-to-noise ratio for detecting and imaging weak components along with strong 
ones. Comparing the data taken with the SMA with the data taken with the eSMA, it was 
found that the inclusion of the CSO and JCMT antennas, i.e. the eSMA, improved the 
signal-to-noise ratio on the longest baselines by 40% compared to that of the SMA alone. 

In order to estimate missing flux of the two HCN transitions in the v=0 and (0, l^'^, 
0) states, we obtained a single-dish spectrum of the two HCN transitions with the JCMT 
on 2008 May 11th (UT). The spectrum was acquired at an elevation near 30°. The main 
beam efficiency was measured to be 65% with a beam size of ~ 19" at 266 GHz, by using 
the continuum emission of Saturn. We used 106.5 K for the main beam temp erature of the 
planet. Note that the deep PH3 absorption feature (jWeisstein fc Serabyn 1994 ) of the planet 
was taken into account when estimating the main beam temperature of the planet at 266 
GHz. Conversion from the scale to flux density scale was achieved with a conversion 
factor of 32 Jy K'^ for the JCMT data. 



3. RESULTS AND DISCUSSION 

3.1. Continuum Component from Circumstellar Dust and the Stellar 

Photosphere 



Continuum data were created from the portions of the spectra which appear to be free 
from line emission. Figure [T] shows visibility amplitude versus uv distance of the continuum 
component originating from the circumstellar dust, and photosphere of the central star, at ~ 
1.1 mm. The data points were averaged over every ~ 5 kA. The falloff of visibility amplitude 
with uv distance demonstrates that the continuum component is spatially extended and 
partially resolved. In Figure 1, one can find that the amplitude is fiattened at 400 kA. This 
means that there are components contributing to the continuum emission which are smaller 
than 27 or 62 AU in size. On the baselines shorter than 400 kA, the amplitude of the 
continuum emission rises as going towards shorter baseline length, which means that there 
are continuum componen ts which we resolved whos e radii > 27 R^. R^ is estimated to be ~ 
500 Rq= 3.5 X 10^3 cm jMen'shchikov et al.l[200ll ). 



It is known that the continuum emission of IRC-l-10216 varies w ith a period of 635 days 
at submillimeter as well as optical wavelengths (IJenness et al.ll2p02n. At 850um. th e average 
fiux is 8.8 Jy and it changes by ± 0.95 Jy due to the pulsation (IJenness et al.ll2002l ). Taking 



^http://www. aips.nrao.edu/ 
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the sinusoidal curve of Figure 7 in iJenness et al.l (|2002| ). at the time when the object was 
observed in 2008 March and April, the object is ex pected to be in the maximu m phase (phase 
of 0.015 by counting from the maximum phase). iGroenewegen et al.l (119971 ) measured that 
the star's continuum emission extends up to 50" at the 3a level of their measurements at 
1.3 mm. The continuum emission measured from our JCMT observation was 3.7 Jy at 1.1 
mm, which i s cons istent with measurements done by other authors (e.g., Dehaes et al. 2007). 



Sahai et al.l (119891 ) estimated the continuum emission from the stellar photosphere would be 



0.55 Jy at 1 mm. It suggests that about half of the total continuum emission imaged with 
the eSMA may be due to the emission of the stellar photosphere, considering the measured 
flux of the continuum component with the SMA was 1.18 Jy (see Section 2). The rest of the 
flux measured with the eSMA may come from circumstellar dust components and low-level 
line emission which were not individually detected. 



3.2. HCN and Other Molecular Species Detected with the eSMA 

Figure [2] shows the spectra of HCN J = 3 — 2 transitions in the v=0 and (0, 1^°, 0) 
states over sky frequency (Fsky) taken at the JCMT and at the eSMA, respectively. Figure 
[3] shows the comparison of the J = 3 — 2, v =(0, l^'', 0) lines observed with the JCMT and 
the eSMA more closely. The eSMA spectra shown in Figure [2] and [3] were produced from 
a spectral cube generated by the image processing. To create the spectrum, a region was 
specifled which includes all emission of the transitions detected above 3a level in the total 
integrated intensity map of the transitions. 



Ziurys fc Turner! (119861 ) observed HCN J = 3 — 2, v={0, 1^'^, 0) transition in 1986 March 
at the Kitt Peak 12m telescope. The stellar phase in 1986 March was near minimum. The 
peak flux observed in 2008 April was ~ 25 % higher than that observed in 1986. The 
increased flux might be due to variability of the stellar luminosity. 

While both the HCN J = 3 — 2, ^;=0 and (0, l^'^, 0) single-dish spectra show a typical 
parabolic-shape proflle that traces the expanding envelope, the line proflles of the transitions 
taken with the interferometer have strange irregular shapes with a few peaky components. 
The fractions of the flux in the f =0 and (0, 1^*^, 0) states detected by the eSMA are only 1 
and 45 % of the emission detected by the JCMT. Please note that the size of the JCMT's 
beam corresponds to 1.1 x 10^ i?* (= 2.6 x 10^ AU). The rest, i.e., the missing components, 
probably originating in the expanding envelope, are resolved out with this array conflgura- 
tion. The size of the missing components is unknown because mapping observation using a 
single-dish telescope has not done. 
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For the line data sets, we discuss only the detected compact components. Physical 
properties estimated from the velocity widths follow in Sections 13.41 and 13. 6[ 

In addition to the two HCN transitions above, seven other molecular transitions are 
clearly detected (Figure 4). One of the detected HCN transitions is in the vibrational state 
of (0, 1^^, 1), mentioned in Section 1. Molecular lines were identified mainly by comparing the 
observed frequencies with the line catalogs provided by NIST (SLAIM), JPL and Universitat 
zu Koln (CDMS). When there were multiple candidates for a spectral line feature, we checked 
the expected line strengths of the different transitions over the observed frequency ranges 
and identified most likely one. For the lines that don't have other transitions in the band, 
we selected most likely one considering expected abundance in the circumstellar region. 

Table [T] summarizes features of nine lines detected with the eSMA. It includes an iso- 
topologue of silicon monosulfide, Si^^S (f=0), potassium chloride, KCl {v = 2), and two 
other HCN transitions in different vibrational states, f=(l, 1^^, 0) and (0, 1^^, 1). The ex- 
citation temperatures of these transitions are about 5800 K and 4000 K, respectively. The 
HCN V = line has a relatively large shift to the red (~ 2.7 MHz or 3 km s~^). This 



may be due to absorption from the hyperfine transitions (IMorris et al.lll985l ) because of the 
nitrogen's quadrapole coupling. The rest of the lines have smaller differences between the 
catalog and observed frequencies (< 700 kHz). The central velocities of molecular transitions 
may be slightly shifted because (1) some of the circumstellar shells are not perfect shells, 
(2) evolution of the molecular species, (3) optical depth of the molecular lines, and so on. 
The excitation of the HCN transitions will be discussed in Section The distributions of 
these hot HCN transitions will be described in Section 13. 4[ 

We detected the J = 35 — 34 transition of KCl in the v = 2 state. KCl was first detected 



at 3 mm by lCernicharo fc GuelinI (119871 ) along with other metals in the circumstellar envelope 
of IRC-l-10216. The transitions they detected were all in the v = state. This is the first 
detection of a KCl transition in a vibrationally excited state. The Si^'^S J = 15 — 14 transition 
has a broader line width (full width at half maximum of ~ 9.2 km s~^) with two peaks. It 
seems to have a wing component on the red-shifted side. The broad line feature of Si^'^S 
might be due to blending from other molecular line(s). However, the lines near the frequency 
were not found in the available catalogs. The line shapes of these transitions will be discussed 
in Section 13. 6[ The maps of these molecular transitions are shown in Section 13.41 

There were four other spectral features whose phases are tightly clustered within ~ 
50° from zero over ^6-10 km s^^, signifying higher likelihood of detection, albeit at much 
lower amplitudes compared to the nine molecular lines described above. Their frequencies of 
265.311, 265.323, 265.347, and 265.449 GHz, are marked as a, b, c and d in Figure H These 
four features are likely detections. 
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3.3. Excitation of the Hot HCN Lines and the KCl Lines 

Excitation of the HCN hnes in three different vibrational states and the KCl line are of 
interest because they exist in the innermost circumstellar envelope where dust and molecular 
species are forming and the molecular gas is being accelerated (see Section [331) . 

The excitation temperatures of the lower energy levels of these transitions, ~ 5800 K 
for the v={l, l", 0) state and ~ 4000 K for the v={0, l", 1) state, are much higher than the 
temperature at the inner edge of the envelope (600 - 1700 K; e.g., Men'shchikov et al. 2001 
and references therein) . It would be difficult to excite the transitions coUisionally unless very 
high density gas (> 10^° cm^^ or so) is present in a large region (> 5"~ 700 AU or ~ 300 
R^,). It's very unlikely that the condition is achieved for this case. The upper level of the 
molecules may be populated directly through the strong radiation of the stellar photosphere. 

Figure E] illustrates energy diagrams of selected transitions of HCN, H2, and CO. The 
observed rotational lines in three different vibrational states are drawn by thin horizontal 
lines and are magnified in the diagram. We propose that some of the HCN lines may be 
pumped through interaction with photons of some transitions of molecular hydrogen, those 
of CO lines, and the stellar radiation. The stellar luminosity of IRC+10216 is ~ 2 x 10^ 
Lq and the peak of the spectral energy distribution (SED) of the object is near 10 /im (e.g. 
Cernicharo et al. 1999). In the near infrared band which is near the peak of the SED, such 
as at 2.5, 3.5, and 14 //m, the molecules may be pumped to the vibrational states of HCN in 
the v={l, l^f, 0), (0, 1^-^.1), (0, ll^ 0), (1, 0°, 0), and (0, 0°, 1) states from the ground level 
to the exited states via the stellar radiation. Figure [5] also shows the excitation temperatures 
of H2 V = 1 state are close to that of the HCN (1, 1^, 0) state and some of the transitions 
have energy levels very similar to that of HCN (1, 1^, 0) state, indicating that a direct and 
rapid photon exchange can occur among the transitions. Furthermore, since CO v = 9 — 7 
transition is at ~ 2.5 /im, such photon exchange can occur with the HCN f = (1, 1^, 0) state 
as well. 

The same excitation method, namely direct pumping from the strong infrared radiation 
at ~14 fim of the stellar photosphere, can explain the excitation of the KCl transition in the 
V = 2 state. 

These hot HCN lines and the KCl line may be maser lines as the stellar radiation must 
excite the transitions. If so, the maser lines of the transitions may not be fully saturated 
because their flux is not high. 
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3.4. Distribution of HCN, SiS, and KCl — the Properties of the Innermost 

Circumstellar Envelope 

Integrated intensity maps of hot HCN J = 3 — 2 emission in the v={l, 1", 0) and w=(0, 
1^^, 1) states are shown in Figure [6] (a) and (b) respectively. The distribution of the emission 
is not spatially resolved. The peak of the emission coincides with the stellar position. The 
radius of the distribution is less than 15 R^,. The compact distribution supports the scenario 
that the pumping arises close to the stellar photosphere or at inner edge of the circumstellar 
envelope. 

Figure [6] (c) and (d) show the integrated intensity maps of the transitions of KCl and 
Si^'^S. Both maps show that the emission originates from a region very close to the central 
star, very similar to the HCN maps of the ^^=(1, 1^^ 0) and v={0, 1^^, 1) states. Both of the 
KCl component and the Si^^S component are not significantly resolved. The Si^'^S emission 
is partially resolved. However, most of the emission originates from a compact region near 
the central star. The compact distribution of the Si^^S J = 15 — 14 v = transition and 
the peaky line features (Figure Hj) may indicate that part of the emission arises from maser 
action. The distribution of KCl v = 2 emission is very compact and the emission is not 
spatially resolved, most likely tracing inner edge of the envelope. 

Total integrated intensity maps of the HCN f=(0, 1^*^, 0) and v = transitions are 
shown in Figure [71 In contrast to the maps shown in Figure [6] that are not spatially resolved, 
the emission is much more extended. Most of the emission is found within a radius of 1", 
which corresponds to ~ 60 R^. One can find two peaks along the northwest - southeast 
direction on both maps. The two peaks of both maps are displaced from the stellar position. 
The position angles {P-A.) of the line drawing between two peak components in the f=(0, 
1^'^, 0) and in the v = states are 115° and 130°, respectively. The distances of the two 
condensations of the ^^=(0, l^'', 0) map are 90 mas and 300 mas from the central star. The 
two condensations of f = map are distributed more or less symmetrically at a distance of 
200 mas from the central star. In the region at a distance of 200 mas from the star, the gas 
temperature is estimated to be a few 100 K (e.g., Keady et al. 1988). 



Murakawa et al.l (120051 ) reported that there are features with a low degree of linear 
polarization identified at 1.5 fxm in continuum emission. The authors concluded that low 
degree of linear polarization found inside the two features indicates that the density in the 
two regions is higher than that of the surrounding region. The two features called NW ellipse 



and SE fan by lMurakawa et al.l (120051 ) may be the same entities as the two HCN peaks traced 
with the HCN J = 3 — 2 transition in the 'i'=(0, 1^'^, 0) state, because of their similarity in 
shapes, P. A., and the fact that both of them trace hig h density regions. However, comparing 



the locations of the polarization minima reported by iMurakawa et al.l (120051 ) and the HCN 
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v={0, 1^*^, 0) peaks in Figure [7] (left) relative to the position of the central star, the HCN 
peaks are located at a distance from the central star only 60% of the distance between the 



polarization minima from the star. Please note that iMurakawa et al.l (120051 ) estimated the 
location of the central star based on their polarization map. It is not clear why the HCN 
peaks and the near infrared polarization minima lie on different locations. If it were due to 
time evolution, the two HCN peaks and the polarization minima might be tracing different 
entities because the HCN peaks are located closer compared to the polarization minima. It 
could be that the HCN peaks detected this time in 2008 might not have existed clearly when 
Murakawa et al made the observation in 2003 January. However, in that case, one can say 
that the direction of the most copious mass loss hasn't changed over about 5 years. 

From the two HCN peaks and those polarization features, we propose that this direction 
along the two HCN peaks of the HCN transitions may correspond to the direction of most 
conspicuous mass loss, and that this direction is the equatorial direct ion of the star. This 



proposition is consistent with the analysis of iMurakawa et al.l (120051 ) who discussed that 
the two linear polarization features may indicate the formation of an edge-on dust torus 
surrounding the carbon star. Discussion on the kinematics of the two HCN peaks and other 
surrounding components will follow in next section. 

Note that there are a few other high angular resolution observations at near infrared 
(e.g., Weigelt et al. 2002; Menut et al. 2007). However, the location of the central star in 
their maps is uncertain so that it is difficult to compare their observed features with ours. 
In particular, we can't find any of likely components that may correspond to the two HCN 
peaks/ the two polarization minima. 



At near infrared (2.2^m), iKastner fc Weintraubl (119941 ) and ISkinner et al.l (119981 ) iden 



tified a 10"-size bipolar reflection nebula that is produced by scattering of radiation from 
the star by dust particles in the region. The axis of the bipolar reflection nebula is at P. A. 
of ~ 30°, which is almost perpendicular to that of the most conspicuous mass- loss direction 
described above. The axis of the bipolar reflection nebula supports the scenario that the 
mass loss is most copious in the direction of the two HCN peaks shown in the maps in Figure 
[7] and that infrared radiation can pass preferentially in the direction perpendicular to the 
most conspicuous mass-loss direction. 

The significant difference in distribution between HCN emission in v= 0/(0, 1^'^, 0) states 
and the HCN emission in the v={l, 1^^, 0)/v={0, 1^^, 1) states or KCl and Si^^S emission leads 
us to conclude that there are two discrete regions in the innermost circumstellar envelope. 
One is a region of radius less than 15 R^^ where various molecular species have just formed and 
the gas is not accelerated by stellar radiation yet. We name this region region I. The other 
is a region outside of the region I, region II, where spectral lines indicate higher expansion 
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velocity. These two discrete regions may have been created because the size distributions 
of the condensed dust particles are significantly different in these regions (see also Section 
13. 6p . Dust nucleation may have further proceeded because temperature decreases rapidly as 
going outwards from the central star. We will discuss the properties of the region II in next 
section. 

3.5. The Kinematics and The Distributions of the HCN v={0, 1^% 0) Maser 

Clumps 

Figure [H] shows channel maps of the f=(0, 1^*^, 0), J = 3 — 2 transition that illustrate 
the detailed internal structures of the gas in region II. Note that the systemic velocity of 
the object is —26 km s~^. The channel maps clearly show expanding shells. The expansion 
velocity of the shells, derived by fitting the spectrum integrated over the whole region, is up 
to ~ 13 km s~^. The expanding shells are composed of many clumpy condensations. 

There are a few prominent condensations along the most copious mass-loss direction. 
By following the maps at different velocities, one can identify shell-like features whose radii 
range from 11 to 22 i?*. Partial shell structures are found in the channel maps at ^sRof 3.8 
and 19.3 km s^^. The velocities are noted in the upper part of each panel. In other words, 
in these velocity ranges, the HCN maser condensations are distributed avoiding the stellar 
position. The radius of the shell is measured to be about ~ 23 /2* with a thickness of ~ 
15 i?*. They may be the smallest shell-like structures identified to date in the circumstellar 
envelope of this star. 

One can find one maser clump in the southeast side of the central star in the channel 
maps of 3.8 - 10.2 km s~^, noted in the upper part of each panel. Another maser clump 
is found in the channel maps of 15.7 - 19.3 km s~^ in the southeast side. Assuming these 
two clumps are located on the east side of the rim of the torus-like feature (see Section 
13. 4p and the clumps have had constant expanding velocities, the medium of the clumps may 
have departed from the central star > 40 years ago. For the northeast side of the central 
star, a clump is seen in the channel maps of 4.7 - 13.8 km s^^, and another one is seen in 
the channel maps of 16.6 - 19.3 km s~^. Assuming a similar condition, i.e., the clumps are 
located on the west side of the rim of the torus-like feature (see Section 13.41) and the clump 
has had constant expanding velocities, the medium of the clumps may have been shed by 
the central star > 100 years ago. 
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3.6. Acceleration of the Gas in the Innermost Circumstellar Envelope 

One can gain insights on the acceleration of the gas due to the stellar radiation from the 
observed line widths and the spatial distribution. The molecular gas is accelerated through 
dust by sharing their momentum through dust-gas collisions as the dust particles receive 
the radiation pressure. Silicon carbide (SiC) can condense into dust particle very close to 
the stellar photosphere at a temperature of ~ 1500 K (e.g. Men'shchikov et al. 2001) and 
they may evolve chemically as going away from the star, while graphite grains form at lower 
temperature of ~ 600 - ~ 1000 K or with radius of 5 - 20 R* (e.g. Men'shchikov et al. 2001). 
At the inner edge of the envelope, SiC may play a major role to accelerate the molecular gas 
outwards. 

One can derive the expansion velocity by fitting line profiles. The line profiles were fitted 
assuming the envelope is spherically symmetric and expanding uniformly with the expansion 
velocity. We also assumed that the turbulence and thermal velocity is much smaller than 
the expanding velocity and the medium is ejected from the central star. Table [T] summarizes 
the expansion velocities and frequency differences between recommended frequencies and 
observed central velocities measured for each line. 

The HCN J = 3 — 2 transitions in the v=0 and (0, 1^°, 0) states have broad line widths 
of ~ 13 km s~^, which is close to the terminal velocity o f ~ 15 km s~^ as observed from 



other molecular lines such as CO (IKnapp fc Morris! Il985l ) and CS (lYoung et al.ll2004l ). In 



region II, both graphite grains and SiC concretion may contribute to the acceleration of the 
molecular gas. The eSMA resolved the emission in both vibrational states. The kinematics 
and distribution of the emission in the (0, l^'', 0) state are investigated in detail (Section 



It is very interesting that Si^^S emission has the largest expansion velocity (~ 16 km 
s~^) although the distribution is compact and is not significantly resolved (in Section [3^ . 
This expansion velocity is at the high end among the expansion velocities summarized in a 
detailed li ne survey study by Patel et al. (2008). Since the sublimation temperature of SiS 



is 1213 K ([Hansen fc Anderkdll962l ). which may correspond to the temperature in an outer 
layer of the stellar photosphere, the SiS concretion mixtured with the SiS gas might receive 
the radiation pressure directly from the stellar radiation, rather than the SiS gas gains the 
momentum from surrounding dust particles in the region. The full width at half maximum 
of the Si^^S line is only 9.2 km s^^. The large expansion velocity is due to the strong wing 
component in the redshifted side. It is known that the sulphur is involved in shock chemistry 
(e.g. Blake et al. 1987). It may arise from a shock region in the innermost envelope. The 
cause of the broad wing emission remains unsolved. 
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On the other hand, the rest of the hnes which have lower expansion velocities of < 7 
km s~^, namely, HCN in the v={l, 1^^, 0) and (1, 1^^, 0) states, the emission of KCl in 
the v=2 state, and all three unidentified lines, must be tracing the very innermost part of 
the envelope because the velocity becomes larger as a funct ion of distance from the stellar 
photosphere (iGoldreich fc Scovilldll976l : IWinters et al.lll994j ). The lowest expansion velocity 
is that for the HCN line in the v={l, 1^^, 0) which has the highest energy level among the 
HCN lines. The small frequency /velocity shifts from the catalog frequencies of the narrow 
line width emission (see Table 1) also supports the scenario that the expansion velocities 
of those components must be very small because their velocities are close to the systemic 
velocity of the star. Even if the molecular components of these transitions have broken shells, 
since the expansion velocities are small, their central velocities should come to very close to 
the systemic velocity of the star. 



4. SUMMARY 

Using the eSMA, we investigated the properties of the innermost circumstellar envelope 
(CSE) of the bright carbon star, IRC+10216. HCN maser emission is spatially resolved for 
the first time in an astronomical object. KCl maser line in a vibrationally excited state is 
detected for the first time. 

We identified two discrete regions in the innermost circumstellar envelope. One is a 
region with a radius of less than 15 i?* where molecular species have just formed and the 
molecular gas has just started being accelerated by stellar radiation (region I). This region is 
traced by hot HCN and KCl lines, and Si^'^S line. Si^^S line detected in region I shows a large 
expansion velocity of 16 km s~^with the broad wing emission unexpectedly. This might be 
partly because SiS is both in gas and solid phase and receives the radiation pressure directly 
from the strong radiation pressure from photosphere. It might be related to some kind 
of shock process in the innermost envelope. The other is a shell region with a radius of 
about 23 and with a thickness of ~ 15 R^, where the expansion velocity has reached up 
to 13 km s^^ (region II). These two discrete regions may have emerged because the size 
distributions of the condensed dust particles are significantly different in these regions, i.e., 
region II may hold larger dust particles or graphite grains. We identified the most copious 
mass loss direction using the HCN J = 3 — 2 transition in the v={0, 1^^, 0) and v = 
states. We propose that this direction corresponds to the equatorial direction of the star. 
The direction was found to be perpendicular to the near infrared bipolar reflection nebula 
and to be parallel to the features with a low degree of linear polarization identified at near 
infrared. The kinematics in region II are studied in detail using the HCN J = 3 — 2 emission 
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in the v—{0, 1^^, 0) state, which shows clumpy structures. Small shell-like features with a 
radius of 11 - 22 i?* were identified in region II. A torus- like feature with a radius of ~9 
exists in region II. Assuming the expansion velocity has been constant, the medium of the 
most conspicuous HCN maser clumps may have been shed by the central star about 40 - 
100 years ago. 

We identified nine molecular transitions in total, including hot HCN J = 3 — 2 maser 
transitions in the ^=(0, 1^^, 1) and v={l, 1^^, 0) states and KCl {v = 2) and Si^^S emission. 
These transitions may be excited directly by the strong infrared radiation from the stellar 
photosphere of the star and by line overlap of some of CO and H2 transitions for the HCN 
transition in the ^^=(1, l^^ 0) state. From the observed line profiles, expansion velocities of 
each line are estimated. All the emission fines except for Si^^S line have expansion velocities 
lower than the terminal velocity of the gas (~ 15 km s~^) in the circumstellar envelope, 
which means that the observed emission traces the region where the acceleration due to 
strong stellar radiation is taking place. The HCN transition in the v={l, 1^^, 0) state has 
the lowest expansion velocity of ~ 3 km s^^, suggesting that it arises from the innermost 
part of the envelope among the lines detected in this observation. Region II has an expansion 
velocity of 13 km s~^. 

Future observations with eSMA will allow us to constrain physical and chemical proper- 
ties of the innermost circumstellar envelope further in detail and construct a concrete picture 
of the mass loss process in a transition phase from AGB to post AGB/proto-planetary nebula 
(PPN). 
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Table 1. Molecular Transitions Detected with the eSMA 









Identification 




t^exp 


Feat — Fcnt 


References for 


(MHz) 


(MHz) 


(Jy) 




(K) 


(km s-l) 


(MHz) 


Frequencies 


U264665 




0.6 






4.2 (6) 






264789.76 


264790.50 


0.8 


Si^^S J = 15 - 14, = 


88.9 


16.3 (2) f 


-0.7 


(1) 


U265253 




0.5 






5.1 (3) 






265264.80 


265265.43 


0.8 


KCl J = 35 - 34, I) = 2 


1012 


4.1 (3) 


- 0.6 


(1) 


265364.36 


265364.90 


1.2 


HCN J = 3 - 2, v={0, l", 1) 


4045 


6.8 (7) 


-0.6 


(2) 


265373.12 


265373.99 


0.6 


HON J = 3 - 2, v={l, l", 0) 


5773 


3.2 (3) 


-0.5 


(2) 


U265386 




0.5 






3.5 (7) 






265852.71 


265852.41 


29 


HCN J = 3- 2, D=(0, l^", 0) 


1036 


12.8 (1) 


+0.3 


(1) 


265886.19 


265883.50 


10 


HCN J = 3 - 2, v=0 


12.8 


12.9 (1) e 


+2.7 


(1) 



^ Recommended frequency from catalog. 
Observed center frequency. 

Peak flux density observed with the interferometer. 

Energy of the lower level of the transition relative to the ground state. 
® Expansion velocity. The numbers in parentheses represent one standard deviation in units of the last signiflcant digit. 
' The line's full width half maximum is 9.2 km s^'^. 
^ The broadening effect due to hyperfine components is subtracted. 

Note. — References: (1) NIST catalog; (2) De Lucia and Helminger 1977. 
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Fig. 1. — Visibility amplitude versus uv distance of the continuum emission of the circum- 
stellar dust and photosphere at 1.1 mm. 
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Fig. 2. — HCN J — 3 — 2 spectrum in the ground vibrational state and ^=(0, l-*^^, 0) state 
taken at the JCMT (top) and with the eSMA (bottom). 
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Fig. 3.— Comparison of HCN J^3-2v =(0, 0) spectra taken with the JCMT (dashed 
hne) and with the eSMA (sohd hne). 
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Fig. 4. — Averaged spectra over all baselines. 
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Fig. 5. — Energy diagram of selected transitions and vibrational states of HCN, H2, and 
CO. The red dots plotted between rotational states marked with red color indicate the 
observed HCN transitions with the eSMA. Under thr ee vibrational mod es (see Section [1]), 
three cartoons are drawn for each vibrational mode (jPichamuthul Il974l ). The blue dotted 
thick arrows pointing upwards and downwards represent excitation and de-excitation of the 
transitions. The orange dotted arrows represent rovibrational maser transitions. The v=9- 
7 transition of CO encased with dotted line at top right is not drawn in the scale of the 
diagram. 
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Fig. 6. — Integrated intensity maps of (a) HCN J = 3 — 2 transitions in the v = (1, 1^^,0) 
state, (b) in the v = (0, l", 1) state, (c) KCl J = 35 - 34 in the ^; = 2 state, (d) Si^^S 
J = 15 — 14 in the v — state overlaid on the continuum emission of the circumstellar dust 
and photosphere. Contours are set at -3, 3, 9, and 15 cr (1 cr ~ 23 mJy). 
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Fig. 7. — Total integrated intensity maps (contours) of HCN J = 3 — 2 v = (0, l^*', 0) (left) 
and V = (right) states, overlaid on the continuum emission of the circumstellar dust and 
photosphere. For the v = (0, l^'', 0) state map, the contours are drawn at -3, 3, 10, 20, 30, 
40, 50, and 60 cr (1 cr = 57 mJy), respectively. For the ground vibrational state map, the 
contours are drawn at -3, 3, 6, and 9 a (1 cr = 39 mJy). 
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Fig. 8. — Channel map (contours) of HCN J = 3 — 2, v=(0, l^*^, 0) transition, overlaid on 
the continuum emission of the circumstellar dust and photosphere. Contours are set at —3, 
3, 9, 15, 27, 33, 39, and 45 a {1 a = 200 mJy). Observed Vlsr is shown at the top of each 
panel. The VisRConverted to that of v={0, 1^°, 0) J = 3 — 2 transition is written in the 
parenthesis at the bottom of each panel. The systemic velocity of the object is —26 km s~^. 



